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Abstract 

Perturbation theory at finite temperature suffers from well-known infrared 
problems. In the standard model, as a result, one cannot calculate the effective 
potential for arbitrarily small values of 0, the Higgs expectation value. Because 
the Higgs field is now known not to be extremely light, it is necessary to determine 
whether perturbation theory is a reliable guide to properties of the weak phase 
transition. In this note, we evaluate the most singular contributions to the po- 
tential at two loops as well as the leading strong interaction contributions. Above 
the critical temperature, the strong interaction corrections are reasonably small, 
while the weak corrections are about 10%, even for rather small values of the Higgs 
field. At the critical temperature, the weak corrections have a more substantial 
effect, rendering the transition significantly more first order, but not significantly 
changing the upper bound on the Higgs mass required for baryogenesis. 
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1. Introduction 



The possibility that the observed baryon asymmetry may have been created at 
the electroweak phase transition has renewed interest in understanding the details 
of this transition. The early studies of Kirzhnitz and Linde [1] and others indicated 
that in the minimal standard model, if m# <C M]y (mji and Myy are the Higgs 
and W mass, respectively) this transition is first order, becoming more weakly 
so as the Higgs mass increases. This result was based on a study of the leading 
terms in the finite temperature potential. More recently, questions were raised 
about this analysis; it was argued that infrared problems might render the phase 
transition more [2] or less [3] weakly first order. In particular, it was argued that 
the potential, at two-loop order, should contain linear terms in the scalar field, 
4> (though of differing signs). More careful study of the perturbation expansion, 
however, shows that these linear terms cancel at this order; the phase transition 
remains first order, though slightly less so [4]. 

Still, the finite temperature theory is full of infrared problems. As a result, it 
is not possible to answer all interesting questions in perturbation theory, even if 
the coupling is weak. The problems are easily illustrated in the minimal standard 
model, with a single Higgs doublet, 0. At high temperatures, the theory becomes 
essentially three-dimensional. For massless particles, the Feynman diagrams of 
such a theory exhibit power-law divergences. In the standard model, the gauge 
boson mass acts as an infrared cutoff. This means that for small 0, the potential 
cannot be reliably calculated. The results of ref. [4] imply that perturbation theory 
is good provided M\y((/)) ^> g 2 T; g here denotes the gauge boson coupling? Had 
the analysis of refs. [2] and [3] been correct, one would have had to impose the 
stronger condition Myj/(0) ^> gT. We are not really interested in the behavior of 
the theory for arbitrary values of the fields, however. One question which we would 
like to address in perturbation theory is the nature of the phase transition. For 
this, one needs to study the behavior of the potential near the minimum. As we will 
review below, for weak coupling in this theory the minimum of the potential occurs 
for g<p oc g T/X, where A is the quartic Higgs coupling. Thus for sufficiently small 
A, the perturbation expansion should be reliable. Indeed, the expansion becomes 
an expansion in powers of g 2 , A, and X/g 2 . If perturbation theory is to be a useful 
guide to the behavior of this system, all of these quantities must be small. 

On the other hand, we already know from LEP that the Higgs particle cannot 
be extremely light, and so the quartic coupling cannot be arbitrarily small. Even 
in extensions of the standard model, designed to produce a more strongly first 



f We will set the Weinberg angle to zero in most of our discussion, since it is not germane to 
the issues at hand; it is easily restored where appropriate. 
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order transition [5] there are likely to be important constraints on the smallness 
of various couplings. Moreover, at finite temperature, one does not necessarily 
obtain the factors of 4n familiar in zero temperature calculations, and corrections 
are potentially large. It is thus natural to ask whether perturbation theory is a 
good guide in the experimentally allowed regions of parameter space. 

In this note, we investigate some of these questions at two-loop order in the 
minimal standard model. It is known that in this model one cannot produce 
an asymmetry, given the present bound on the Higgs mass [4,6]. We choose the 
minimal model to illustrate our points because of its simplicity; our analysis is 
readily extended (often trivially extended) to other theories in which the phase 
transition is likely to be more strongly first order. A complete evaluation of all 
two-loop diagrams is quite involved. However, it turns out that there are certain 
leading contributions which can be easily computed. The simplest to describe are 
the QCD corrections. These are potentially important for top quark loops, due to 
the large value of the top quark mass and of the strong gauge coupling. They can 
be obtained by suitably adapting existing QED computations. In addition, there 
are certain purely electroweak effects which are enhanced by powers of logarithms 
of the scalar field. We have already remarked that in ref. [4] it was shown that 
certain leading infrared effects cancel. It turns out that the subleading terms come 
with logarithms of the gauge boson mass over the temperature. We will compute 
a set of terms in the potential proportional to g 4 2 Iog(g0/2T). We will see that, 
away from the critical temperature, T c , for the phase transition, both the QCD and 
pure electroweak effects give corrections of order a few per cent to the potential. 
Thus perturbation theory appears to be in quite good shape. Near T c , however, 
the situation is more delicate; the weak interaction corrections are in fact quite 
large, here, and seem to make the transition more first order. 

2. Review of Earlier Results 

The zero temperature potential, taking into account one-loop corrections, is 
given by [7,8] 

V = "^0 2 + \^ + 2Br$<? - + In (£) . (2.1) 

Here 

B = 64^ ( 2m W +™ 4 z- Ami) , (2-2) 
v = 246 GeV is the value of the scalar field at the minimum of Vq, A = jJ? /v 2 , 
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m\ = 2/i 2 . At a finite temperature, one should add to this expression the term 

v T = £2 W-{vw) + u-(yz) - ui+ivt)] , (2.3) 

where m = Mi<p/v T, and 

00 

I T (y) = ± J dx x 2 In T e-V^^j . (2.4) 


In the high temperature limit it is sufficient to use an approximate expression for 
V(<j>,T) [1,5] 

V(<j>, T) = D{T 2 - T 2 )(j) 2 - ETcj) 3 + . (2.5) 

Here 

D = ± (2m 2 w + m 2 z + 2m 2 ) , (2.6) 
E = -L (2vr? w + m%) ~ 1(T 2 , (2.7) 

T ° = 2^ (/i2 " ABV ° ] = 4^ (m ^ ~ 8BV ° ] ' (2 ' 8) 

AT = A-^f 2 m^ln^ + m|ln^-4m|ln^|,) , (2.9) 

where lnas = 21n47r — 27 ~ 3.91, lnap = 21n7r — 27 ~ 1.14. As explained in 
refs.[4] and [9], the constant E is reduced by higher order corrections by a factor 
of 2/3. 

This potential leads to a phase transition which is at least weakly first order, 
basically as a consequence of the term cubic in 0. At very high temperatures, 
V has a unique minimum at = 0. As the temperature is decreased, a second 
minimum appears at a temperature 

T * = i-Jw - (2 ' 10) 

At a temperature T c , this new minimum becomes degenerate with the minimum 
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Fig. 1. Two- loop contribution to the potential involving a gluon correction to a top 
quark loop. 



at the origin; T c and the corresponding minimum of the potential C are given by 

T 2 - S fk - 2ETc (oil) 

Finally, at T , the minimum at the origin disappears, and the potential has a 
unique minimum. 



3. QCD Corrections 



Because they are conceptually simplest, we consider, first, the strong correc- 
tions to the potential. The leading terms arise from a gluonic correction to a top 
quark loop (fig. 1.) Apart from a trivial group theory factor, this calculation is 
identical to the calculation of the two- loop free energy in QED [13,14]. 

If one examines the Feynman diagram containing the top quark (electron in 
QED) they are completely free of infrared problems even as the quark (electron) 
mass tends to zero. This is because the discrete Matsubara frequencies for fermions 
are all non-zero. However, the results of the QED calculation [13,14] diverge loga- 
rithmically with the mass [as m 2 , ln(m e )]. The reason is not difficult to understand. 
The QED calculation is performed by subtracting the electron mass on shell. The 
required counterterm itself diverges logarithmically. In QCD, we want to subtract 
the mass at some off-shell point; in the present case, we are interested in the mass 
(or, more precisely, gt), at a scale of order T, the temperature. Renormalizing off 
mass shell in this way, there are no logarithms. 

To evaluate the top quark contribution, we begin with the QED expression (for 
zero chemical potential) [13,14] 



V 



; 2 T 2 



d 3 p n F (p) 
(2tt) 3 E p 



d 3 p d 3 q 1 
2tt 3 2tt 3 EpE q 



x 2 
where 



1 + 



m 



+ 



m 



(E p -E q y-{p-q)* ( Ep + E q )2-(p-q)2 

1 



77i 2 * 2 i 2 

E = p + m 



n F (p) 



n F (p)n F (q) (3.1) 



e l3E v + i ■ 

We will content ourselves with extracting the term of order m 2 . The terms 
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involving involving m 2 explicitly in eq. (3.1) must be handled with some care (we 
thank Peter Arnold for pointing this out to us). It is necessary to perform the 
angular integrals first before taking the m — > limit. The integration is then 
straightforward, apart from the integral 



/ 



dxdy ln [ X y \l - i = — (ln(2) - 1) (3.2) 

y (x + y) 2 (1 + e x )(l +ev) 3 K K ' ' K J 



(We thank Howard Haber for evaluation of this integral.) The remaining integrals 
factorize, and are readily evaluated using standard tricks (e.g., , see appendix A 
of ref. [14]). In particular, the integral 

oo 

1 f x 2 dx 1 

= o / t 2 i 2U/2 ( 3 - 3 ) 

2 J (x 1 + y z yl z P \/x 2 +y 2 , 1 



is given by 



_2 



/3 = ^ + y (lny + 7-lnvr-i) . (3.4) 



As we have noted above, this does not yet give us the result we want; we need to 
add a finite counterterm to change the result from an on-mass-shell subtraction to 
an off-shell one. We choose to use the MS scheme. A straightforward calculation 
gives that in the one-loop result, one must replace m by 



m — > m < 1 — 



3e 



2 



In 



m 2 \ 4 



(3.5) 



To order m 2 , the lowest order result is m 2 /T 2 . Putting everything together then 
gives for the fermionic contribution 

-7tt 2 T 4 m 2 T 2 5g 2 T 4 g 2 m 2 T 2 r , , . . 2l , ft . n . . 

^/ = — ^- + ^ + ^- + ^2-h|- ln W+7 + |ln(2)] . (3.6) 

The corrections both to the leading T 4 term and to the quadratic term are about 
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Fig. 2. Two-loop diagrams which potentially give linear terms in the potential. 



Fig. 3. The two diagrams of fig. 2 can be viewed as a correction to the one-loop gauge 
boson propagator. 

4. Pure Electroweak Corrections 

We now turn to the purely electroweak corrections. It is important, first, to 
understand the nature of the perturbation expansion at high temperatures. In 
general, if one examines the Feynman diagrams for the potential, one sees that 
they diverge in the infrared as the gauge boson masses (~ g<f>/2) tend to zero. The 
divergences are cut off by the gauge boson masses themselves. For example, at 
one-loop, the coefficient of the quartic term in the potential is linearly divergent 
as the mass tends to zero; this is the origin of the cubic term in the potential. At 
higher loops one encounters more and more severe divergences. From the start, 
then, it is clear that perturbation theory cannot be good for arbitrarily small values 
of 0. One can hope, however, that if A is small enough, the expectation value of 
the scalar field will be large enough that perturbation theory will be valid near the 
minimum of the potential. Note that the one-loop potential gives for the location 
of the minimum at T c a result which goes as g 3 / A. The potential at the minimum 
is of order g 12 /\ 3 . 

The analysis of refs. [2] and [3] suggested that at two-loop order the coefficient 
of the quadratic term in the potential is linearly divergent. Extended to higher 
orders, this analysis would suggest that the loop expansion parameter is g 2 /My^; 
for oc g 3 /A, this gives \ 2 /g 6 at the minimum. In ref. [4], it was shown that there 
are no such linear terms; a simple extension of the arguments given there shows 
that the expansion parameter is g 2 /M\y, or \/g 2 at the minimum. 

Actually, while it was shown in ref. [4] that the leading divergences cancel, the 
subleading terms come with logarithms of the gauge boson mass, i.e., with ln(0 2 ). 
To understand how this comes about, it is helpful to describe the results of ref. [4] in 
a slightly different fashion than explained in that paper. We assume that the Higgs 
mass is small and neglect Higgs self-interactions. We work, as there, in Coulomb 
gauge; the theory in the infrared is then like a three-dimensional theory in Landau 
gauge. As explained there, for small one can ignore the contributions of the 
longitudinal modes. The diagrams involving transverse gauge boson loops of fig. 2 
are each potentially quite singular in the infrared. Indeed, while each diagram is 
formally of order 2 , power-counting suggests that the integrals separately behave 
as 1/0, and that these diagrams should thus give linear terms in (this is the 
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source of the linear terms in refs. [2] and [3]). However, the leading divergences 
cancel between the two diagrams. This is easily understood; the diagrams can be 
organized as in fig. 3. In other words, one can view the combination of diagrams as 
a propagator correction to the one-loop gauge boson contribution to the free energy. 
To verify this requires care with the combinatorics; as explained in ref. [4], it is 
simplest to verify statements of this kind for tadpoles, rather than directly for the 
free energy. In any case, it is well-known that the three-dimensional gauge boson 
propagator vanishes at zero momentum as a consequence of gauge invariance. A 
straightforward calculation gives 

32 a/ q z 

Because the polarization tensor vanishes linearly with the momentum, the integral 
behaves for small as ln(0). A simple calculation gives 

* l/l = - 2 li? lnW/2r) - (4 - 2) 

(note that, as is conventional in the finite temperature theory, we are again nor- 
malizing clS du real field here). 

There are other corrections at two loop order which depend on log(0). These 
come from the diagrams in fig. 4, which contain scalar fields in intermediate states. 
These are somewhat more subtle to treat, since one must perform a resummation of 
some type in order to determine the effective scalar propagator. However, without 
probing this issue too deeply, it seems likely that any resummed scalar mass will 
be smaller than the gauge boson mass. For example, even well above the phase 
transition, the effective scalar mass in the one loop potential is somewhat smaller 
than the gauge boson mass; near the transition, it is much smaller. So to get an 
estimate of these effects, we will simply neglect it. In this case, the diagrams are 
straightforward to evaluate? From the first, we obtain 

SV 2 = ^^H g( p/2T), (4.3) 



★ This calculation can be found in many textbooks; see, e.g., , refs. [16-17]. However, note 
that the formula the first reference contains a typographical error; we have verified that the 
formula in the second is correct. 

f In evaluating the first diagram, one approach is to route the loop momenta, I and k so 
that the scalar propagator depends on k — I. This propagator can then be expanded in a 
(double) series of Legendre polynomials. The momentum integrals become trivial, and the 
sums give (-functions. 
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Fig. 4. Diagrams with scalers in intermediate states. 



while from the second we have 



SV 3 = 



35W 
256tt 2 



ln(#/2T) . 



(4.4) 



So overall, the shift is 



6V1 = 



51# 4 T 2 2 

5127T 2 



ln(#/2T) . 



(4.5) 



Away from the critical temperature, this correction is not too large. If the 
logarithm is of order 2 (corresponding to typical 0's at the transition), we only 
obtain a 10% correction. However, near the phase transition, the effect of this 
term is enhanced, because the quadratic term is so small. For example, keeping 
only the one-loop correction, the coefficient of the quadratic term at T c is of order 
10~ 2 T 2 . This is smaller than the size of the logarithmic term. So the effect of the 
logarithmic correction is potentially dramatic. In order to evaluate its significance, 
note first that the various two-loop corrections we have considered alter the value 
of the critical temperature. Taking mjj = 35 GeV, we have first found this new 
temperature, T c In fig. 5, we have presented two curves: the one-loop result (with 
the correction of ref. [4]) at T c , and the two-loop result, at T c . The two-loop cor- 
rection is seen to render the transition more first order; at the peak, the correction 
is almost a factor of 3. One can see in fig. 5 that the location of the minimum 
changes by about 20%. This weakens somewhat (by about 4 or 5 GeV) the bounds 
mentioned earlier on the Higgs mass. 

It is reasonable to ask about the gauge-dependence of this result. This is a 
subtle question which we will not explore here. We expect, however, that for 
physically meaningful questions (such as the sphaleron rate at the minimum) these 
results should be gauge-independent. 
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Fig. 5. One-loop potential at T c and two- loop potential at T c for a 35 GeV Higgs. 

5. Conclusions 

From this work, it appears that perturbation theory at finite temperature in 
the standard model is not in terribly bad shape, even for Higgs fields which are 
not extremely light. Away from the critical temperature, perturbative corrections 
are quite small. Near T c , however, the situation is somewhat more subtle. Because 
the leading quadratic terms in the potential vanish at this point, subleading terms 
which are not precisely quadratic (i.e., have logarithmic modifications) are of great 
importance. At two-loop order, in the minimal standard model, they render the 
transition significantly more first order, but only slightly modify the limits on the 
Higgs mass. 

It is natural to ask about the effects of still higher order corrections. Examina- 
tion of some particular graphs leads us to suspect that the picture presented here 
will not be drastically altered. Still, given the size of the corrections which we have 
found, non-perturbative calculations would be of great interest. In particular, we 
have not carefully investigated the problem of corrections to the potential loops 
with Higgs fields only, which might be quite important at the transition. 
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